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Abstract 

Based on a rigorous quantum-kinetic approach, spin-charge coupled drift-diffusion equations are 
derived for a strongly confined two-dimensional hole gas. An electric field leads to a coupling 
between the spin and charge degrees of freedom. For weak spin-orbit interaction, this coupling 
gives rise to the intrinsic spin-Ball effect. There exists a threshold value of the spin-orbit coupling 
constant that separates spin diffusion from ballistic spin transport. In the latter regime, undamped 
spin-coherent oscillations are observed. This result is confirmed by an exact microscopic approach 
valid in the ballistic regime. 



PACS numbers: 72.25.-b, 72.10.-d, 72.15.Gd 



1 effect (SHE) (J Q has recently 
reveal an electric-field induced 



The generation and manipulation of a spin polarization by exclusively electronic means 
in nonmagnetic semiconductors at room temperature is a major challenge oi spintronics. 
Among many interesting phenomena, the intrinsic spin-Hal 
attracted considerable interest. Experimental studies 0, 4], 
spin accumulation near the edges of a confined two-dimensional electron (hole) gas. Most 
theoretical interpretations of these experimental data rely on the notion of a spin current 
oriented transverse to the applied electric field. Interestingly, this seemingly clear 

physical picture still remains the subject of serious debates, jg, 7\ The relationship between 
the spin current and the induced spin polarization seems to be a very subtle issue. The 
main problem underlying the debates is the notion of a spin current itself because spin 
is not a conserved quantity in spin-orbit coupled systems. Consequently, any approach 
that avoids the intricate identification of a more or less suitable spin current is superior. 
Such an alternative approach not only introduces different calculational techniques, but also 
suggests alternative interpretations of the effects under consideration. As widely anticipated, 
a complete physical description of spin-related phenomena is provided by microscopic models 
based on the spin-density matrix or Keldysh Green functions together with an analysis of its 
long- wavelength and low-frequency limit. These approaches are more general and free from 
artefact associated with ambiguous definitions of the spin current. In addition, spin-charge 
coupled kinetic equations allow the treatment of such interesting phenomena as propagating 
spin excitations or the relationship between the intrinsic SHE and the zitterbewegung. 

In this report, we propose an alternative approach to the SHE by deriving spin-charge 
coupled drift-diffusion equations for a two-dimensional hole gas (2DHG), which refers to the 
populated heavy-hole band of thin p-type quantum wells. The related heavy-hole Hamilto- 
nian of the cubic Rashba model has the second quantized form 

H = a\ x [e k - e F ] a kX - 2J (^fe * °aa') «L a fcA' 

fe,A fc,A,A' 

+ U ^2Y1 a k\ a k'\ ~ ™E V ' na{_ K/2X a k+tl/ 2x U=o, (1) 

k,k> A fc,A 

where a kX (g^a) denote the creation (annihilation) operators with in-plane quasi-momentum 
k = (k x , k y , 0) and spin A. The electric field E is oriented along the x axis. Furthermore, 
Ef denotes the Fermi energy, a the vector of Pauli matrices, e k = h 2 k 2 /(2m), and u the 
strength of the 'white-noise' elastic impurity scattering, which gives rise to the momentum 
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relaxation time r. Contrary to a phenomenological approach, we treat elastic scattering on 
a full microscopic scope. The spin-orbit coupling is given by 

nw h ^^[i(k 3 + -k 3 _),(k 3 + + k 3 _),o], (2) 

where k± = k x ± ik y and fiuJk = ak 3 . Based on the Born approximation for elastic impurity 
scattering, the Laplace-transformed kinetic equations for the physical components of the 
spin-density matrix have the form [?| 

sf--( K - k)f + i W(S (fc) ■ f + ^V k f = -(f-f) + f , (3) 
m a t 

sf + 2K x /)--(«• fc)/ + iu> K {k)f + ^V fc / 
m n 

where an additional frequency appears 

bJn{k) = —7- \_{ky k x )Ky 2k x kyK x , 

(k x — k y )n x — 2k x kyK y , 0] , (5) 

which depends on k. The cross line over fc-dependent quantities denotes an integration 
over the polar angle ip of the in-plane vector k. f and fo denote the initial charge and spin 
density components, respectively. The quantum Boltzmann equations (j3J) and (jlj) are treated 
in the long-wavelength limit in order to derive spin-charge coupled drift-diffusion equations. 
To this end, the kinetic equations are written in a matrix form Af + Sf = Bf + 5, where 
the matrix A collects all contributions that are independent of the electric field E and not 
integrated over the angle ip [f denotes the four component vector (/, /)]. To calculate the 
matrix B on the right-hand side of this equation, we assume akp -C e^p and restrict k 
contributions up to k 2 . The matrix equation is solved iteratively in the case of weak electric 
field contributions 8 (the matrix 8 contains first-order derivatives and / is decomposed 
according to f + fx with f ~ E° and fx ~ E) . The solution of the equation is written in 
the form 

(l + C- 1 1 C 2 )- 1 C- 1 1 J=l (6) 

where f = C\5 and / 1 = C2S. The general expressions for the transport coefficients are 
very cumbersome but simplify considerably in the low-field case and under the condition 
akp < e kp . 



As we are mainly interested in the SHE, let us focus on the coupling between the spin and 
charge degrees of freedom. By applying the outlined schema, we obtain our main analytical 
result 

(s + i—K x + D k 2 )J + iT z K y J z = /o, (7) 
°"o 

(s + — + i—k x + D z k 2 )] z + iToKyf = 0, (8) 

Tsz &0 

with the transport coefficients 

9Q 2 3a 2 -AQ 2 1 AQ 2 

T °- Vd ^(a 2 +^r (10) 

9Q 2 2 4a tt 2 + 8fi 2 - 3o$st 
z ~ Vd ^ a ° (a sr + 4W)(a 2 Q + 40 2 ) ' ( } 

The parameters introduced in these equations are given by: v<i = eEr/m, <jq = st + 1, 
D = v 2 t/2, 7 = EkT/h, and Q = lo^t = K(k)l (with K(k) = amk 2 /h 2 and the mean- 
free path I = vt). As the contributions ~ k x do not affect our analysis, we took them to 
lowest order in Q. The Eqs. (j7|) and (jSj) have been derived for n/7 < 1 but unrestricted 
values Q. In the absence of the electric field (E = 0), the Eqs. ([7]) and (jHJ) completely 
decouple. This decoupling, which applies to all components of the spin-density matrix, 
represents a speciality of the cubic Rashba model. 0, [^j The situation is completely different 
for electrons (linear Rashba model), for which at zero external fields the charge density / 
couples to the transverse spin component f r = i(K,xf) z , whereas f z is connected with the 
longitudinal component f d = in ■ f. However, both for the linear and cubic Rashba model 
additional couplings arise due to an applied electric field. The related magnetization gives 
rise to a magnetoelectric effect that has been thoroughly investigated in the literature for 
semiconductors with spin-orbit interaction. 

The time dependence of the coupled spin-charge transport is calculated by an inverse 
Laplace transformation of the solution of Eqs. (171) and (El) • Due to the complicated s depen- 
dence of all transport coefficients, a non-Markovian temporal evolution is expected. How- 
ever, the straightforward determination of this complicated time dependence of charge and 
spin densities becomes problematic. Eqs. (J7J) and (|HJ) are only solvable by inverse Fourier- 
and Laplace transformations under appropriate additional approximations concerning the 
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s dependence. This delicate mathematical problem will be accounted for in a forthcoming 
paper. Here, we focus on steady-state solutions (<7o = 1). 

The character of the coupled spin-charge transport strongly depends on the strength of 
the spin-orbit coupling, which is expressed by the dimensionless parameter Q = K{k)l. It is 
the most striking feature of the drift-diffusion Eqs. ([7]) and ([H]) that the character of the spin 
transport changes radically with increasing coupling strength Q. The appearance of such a 
crossover is due to the unusual expression for the diffusion coefficient D z in Eq. (Q, which 
has recently been derived by an alternative approach. [9] Moreover, the very same result is 
also obtained for the linear Rashba model, when the frequency uo^ is appropriately redefined. 
With increasing spin-orbit coupling K(k) or relaxation time r, the diffusion coefficient D z 
changes its sign. A negative diffusion coefficient indicates an instability of the spin system. 
Under this condition, spin diffusion has the tendency to strengthen initial spin fluctuations. 
The competition between this self strengthening and spin relaxation processes results in a 
spatial oscillatory spin pattern. Going from weak (Q < 1/ y/l2) to strong (Q > 1/ VT2) 
spin-orbit coupling, we observe a transition in the spin system from a diffusive behavior to 
a ballistic regime. We shall show that ballistic spin transport is characterized by undamped 
spin oscillations. 

To be specific, let us solve Eqs. ((Tj) and (|HJ) for a stripe geometry (— L < y < L ). The 
inverse Fourier transformation is accomplished by the replacement n y — > id/dy, whereas 
along the x axis all quantities are constant. The resulting set of differential equations is 
easily solved. For the boundary condition /(— L ) = f(L ) = f and f z (—L ) = f z (L ) = 0, 
we obtain 

~f i ; \ , T SZ V ysinh(L // ) - L smh(y/l ) 
D sinh(L //o) 
where / = y/D z r sz . An interesting effect, which we do not follow up in this paper, results 
from the y dependence of the charge density / that is strongly affected by the boundary 
condition and that gives rise to a self-consistent electric field oriented along the y direction. 
In Fig. 1, the thick line illustrates the result for the spin polarization in the diffusive regime 
{D z > 0), when the spin-orbit coupling is weak f2 < l/^/l2. The electric field aligned along 
the x axis induces a spin polarization at the edges of the stripe. This SHE has received a 
great deal of attention. Most popular is the description by means of a spin current oriented 
perpendicular to the electric field. Many theoretical studies (see, for instance, Ref. [9! [lo|) 
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FIG. 1: Out-of plane spin polarization induced by an electric field applied parallel to the stripe 
of a 2DHG (with 7 = 0.5). The thick and thin lines were calculated with f2 = 0.25 and f2 = 0.5, 
respectively. 

introduced the spin current by a symmetrized product of spin and velocity operators 
(v fJ ,a t + a' t v ll )/2. It was claimed that at least for the cubic Rashba model, the SHE introduced 
in such a manner is robust against disorder. Experimental results seem to confirm this 
physical picture. However, there is a principal difficulty with such an approach. The above 
mention definition of the current is not sufficiently general. It would completely fail for any 
hopping transport problem, for which the eigenstates have no dispersion. This definition 
only applies, whenever the Hamiltonian commutes with the dipole operator. Obviously, this 
is not the case for the Rashba Hamiltonian. Consequently, it is necessary to go back to 
the more general definition, which expresses the current by the time derivative of the dipole 
operator. Q] This more general treatment of the spin current reveals a close relationship 
between the field-induced spin accumulation and the spin current expressed by a quasi- 
chemical potential. From a technical point of view, the current that applies in its most 
general form to a homogeneous system is calculated from the quantity V K /(fe, k) at k, = 
and not only from the density matrix f(k). Based on this general definition of the spin 
current, it was concluded that an electric-field-induced steady-state spin-Hall current does 



111 ] On the other hand, a SHE was demonstrated by 



not exist in the cubic Rashba model, 
a recent experiment. This calamity indicates that the notion of a spin current is not 
useful for studying the SHE. An alternative, which is proposed in this paper, determines the 
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spin accumulation from quantum kinetic equations or the associated spin-charge coupled 
drift-diffusion equations. The specific difficulty of the latter approach is the formulation of 
appropriate boundary conditions. Q] 

The character of the SHE dramatically changes in strongly spin-orbit coupled systems 
(Q > I/a/12 so that D z < 0). The result is illustrated in Fig. 1 by the thin line. A spin- 
coherent standing wave travels through the stripe. It is remarkable that these oscillations 
are not damped although a finite elastic scattering is present. The occurrence of a per iodic 



spin pattern is not unusual and has been investigated in the literature. [14|, [15|, [16|, [17( The 
novelty here is that such an oscillatory spin pattern can be induced by an electric field. The 
rapid variation of the out-of plane spin polarization induces a magnetic field that leads to 
circulating microscopic currents. The retroaction of these currents on spin may result in a 
finite damping of spin oscillations. 

In the strong-coupling limit (Q = K(k)l ^> 1), the oscillatory spin pattern changes on 
a length scale K~ 1 (k) that is much smaller than the mean-free path /. This fact conflicts 
with basic assumptions of the drift-diffusion approach, which is only applicable for diffusion 
lengths much smaller than the mean-free path. Although macroscopic transport equations 
were found to be valid even when the spin-diffusion length is somewhat less than the mean- 



free path [18], it is indispensable to treat this point in detail. Large values for f2 give rise 
to spin-relaxation times t sz , which are much smaller than the elastic scattering time r. 
This condition characterizes the ballistic spin regime. Therefore, we go back to the kinetic 
Eqs. ([3]) and (HI) and solve them under the condition r — > oo and to first order in the electric 
field. For the out-of plane spin polarization, we obtain 

(a 2 + 4:LU 2 k )f z + e -^~-^-f z = ^(UkxU K y - UkyUKx)h, (13) 

where a = s — in ■ v. Calculating the inverse Laplace and Fourier transformation and 
integrating over the angle (p, we arrive at the analytical solution 



fz(k,r,t) 



eE fsm2K{k)r 

' cos 2K{k)r 



167rr dk \_Sk \ 2K(k)r 
5(y - a/ (vt) 2 - x 2 ) for 0<y <vt 
-5(y + a/ (vt) 2 - x 2 ) for vt < y < 

which describes the field-induced spin polarization that occurs, when initially a drop of 
charge carriers is injected into the 2DHG at the position r = 0. A 5-like wave front of 



FIG. 2: The amplitude F(r) = r/xd((sin x/x — cos x)/x)/dx in dependence on r with x = 1K(kp)r 
and K(kp) = 0.5. For the factor 77, we have: 77 = ±1 for y ^ 0. 



spin polarization travels through the homogeneous 2DHG. The amplitude of this narrow 
wave front oscillates as illustrated in Fig. 2. As expected, the relief is antisymmetric with 
respect to the y axis. Most interesting for our comparison with the above drift-diffusion 
approach is the observation that the wavelength of the spatial spin pattern is of the order 
of K~ 1 (kp). Therefore, this exact analytical result confirms the existence of a field-induced 
long-living spin pattern in strongly spin-orbit coupled systems as predicted by the drift- 
diffusion equations. 

Our study of electric-field induced spin phenomena revealed a close relationship between 
the SHE and spin-coherent oscillations. We compare this conclusion with recent results that 
demonstrated that the intrinsic SHE and the zitterbewegung are essentially the same kind of 
phenomena. [l?J] Consequently, the question arises whether the above treated spin-coherent 
waves have to be identified with the zitterbewegung, which is a purely relativistic effect. In- 
deed, both kinds of oscillatory spin excitations exhibit common features. The characteristic 
wavelength of both types of oscillations amounts about 100 nm (calculated by adop ting the 
typical values am/h 2 ~ 2 nm and kp ~ 0.1 nm). Moreover, both the zitterbewegung 19| and 
the spin-coherent waves [cf. Eq. ( 1T2|) ] are resonantly enhanced, whenever the width of the 
stripe matches a characteristic wavelength of the spin excitation. However, some features of 
spin-coherent waves are not compatible with such an identification with the zitterbewegung. 
The dispersion relation of coupled spin-charge excitations is calculated from the vanishing 
determinant of Eqs. ([7]) and ([8]). In general, one obtains not only spin-coherent solutions 
but also damped excitations, whereas the transition between them could be driven by the 
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electric field. In addition, the spin-coherent waves that appear at strong spin-orbit coupling 
are separated from the SHE by a sharp threshold. For the zitterbewegung such a threshold 
is not expected as its existence is solely due to at least two energy bands separated by a 
nonzero gap. We think that the interesting study of the relationship between spin-coherent 
waves, the zitterbewegung, and long-living spin-coherent states [fjj] will continue in the near 
future. 

The experimental observation of the field-induced spin-coherent waves should be possible 
by high-resolution scanning-probe microscopy imaging techniques. The direct experimental 
proof of this effect would fascilitate developments both in spintronics and basic research. 

This work was supported by the Deutsche Forschungsgemeinschaft and the Russian Foun- 
dation of Basic Research under the grant number 05-02-04004. 
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